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We study the breakup of a liquid jet a few nanometers in 
diameter, based on a stochastic differential equation derived 
recently by Moseler and Landman [Science 289, 1165 (2000)]. 
In agreement with their simulations, we confirm that noise 
qualitatively changes the characteristics of breakup, leading 
to symmetric profiles. Using the path integral description, 
we find a self-similar profile that describes the most probable 
breakup mode. As illustrated by a simple physical argument, 
noise is the driving force behind pinching, speeding up the 
breakup to make surface tension irrelevant. 



Given the current trend for miniaturization, it is nat- 
ural to ask about the behavior of liquids on increasingly 
small scales, where thermal fluctuations eventually be- 
come important. In the case of free surfaces, relevant for 
all printing and microstructuring applications 0j , fluctu- 
ations can expected to be even larger than in confined 
geometries. Their importance is estimated by comparing 
the thermal energy ksT with the surface tension 7. Thus 
for structures whose size is in the order of the thermal 
length scale £t = {ksT '/^/) 1 ^ 2 0, usually about a nm at 
room temperature, fluctuations should make a leading- 
order contribution. Predicting the (typical) behavior of 
fluids on the nanoscale is therefore an extremely chal- 
lenging statistical physics problem, since the noise can 
no longer be treated in a standard perturbative fashion 
||. Instead, non-perturbative methods are called for to 
properly account for thermal fluctuations. 

In the absence of experiments that could directly mea- 
sure fluid flow on the nanoscale, molecular dynamics 
(MD) [QJ|] computations are an ideal tool. Previous anal- 
yses of drops and threads convincingly showed hy- 
drodynamic behavior on the nanometer scale, and agree- 
ment with breakup times and growth rates as predicted 
by linear theory ||. Recently pJ] a realistic molecu- 
lar dynamics simulation of a jet of propane issuing from 
a nozzle 6 nm in diameter was performed, which also 
payed close attention to the nonlinear dynamics close 
to breakup. Remarkably, a coherent jet develops, from 
which drops separate at slightly irregular intervals to 
form a quasi-stationary decaying state. If nozzles of this 
size were to be built, this opens dazzling perspectives of 
transporting current printing and structuring techniques 
into the nanoworld. 

However, a careful analysis of the breakup process re- 
vealed Jl(J that the presence of noise qualitatively alters 
the character of the breakup. While the deterministic 



motion at a corresponding Reynolds number forms elon- 
gated necks between successive drops M, noise leads to 
profiles symmetric around the pinch point. Thus satel- 
lite drops rarely form and quite surprisingly the distribu- 
tion of drop sizes becomes narrower. In addition, on the 
nanoscale the motion of the minimum neck radius accel- 
erates as breakup is approached, while the correspond- 
ing time dependence h m %n — 0.03(7p/^)(io — t) is linear 
for deterministic pinching [ [TT|| . Here v is the kinematic 
viscosity and p the density of the fluid. Thus the theo- 
retical challenge is to understand this qualitative change 
of behavior in a regime where noise makes the leading 
contribution. 

To deal with the above set of problems, we use the con- 
tinuum description given by ]l0[ | , which consists in adding 
a stochastic term to the lubrication description of a de- 



terministic jet |12|. The amplitude of the noise is fixed 
by the condition of detailed balance Detailed nu- 

merical simulations of the stochastic equation gave very 
convincing agreement with MD simulations. This means 
that hydrodynamics, at least when suitably generalized 
to include fluctuations, is fully capable of describing free- 
surface flows down to a scale of nanometers. 

The coupled set of differential equations for the radius 
of the fluid jet h(z,t) and the axial velocity v(z,t), as 
derived in |1CJ], reads 



d t h 2 + {h 2 v)' = 



(1) 



d t (h 2 v) + (h 2 v 2 )' = -G' + 3(h 2 vJ + D(hg)', (2) 

where the prime refers to differentiation with respect to 
the spatial variable. The first equation (||) expresses mass 
conservation, (^) is the momentum balance. All quanti- 
ties are written in units of the intrinsic length and time 
scales i v — v 2 p/"f, t v = v 3 p 2 /-f 2 , respectively. For later 
convenience the Laplace pressure term is written in the 
form G' — h 2 n' , k being the mean curvature of the in- 
terface. The Gaussian Langevin force is uncorrelated in 
space and time, 



< C(*l,*l)£(<S2,*2) >= 5(zi - z 2 )8(t 1 -t 2 ), 



(3) 



and the dimensionless constant D = (6/tt)£t/(v mea- 
sures the strength of the noise contribution. 

Since the derivative of the noise is an extremely sin- 
gular quantity, it is useful to integrate @ once, setting 
P = f Q p(x)dx, where p = h 2 v is the momentum. Thus 
we arrive at the conserved form of the equations 



d t h 2 



-P" 



(4) 
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d t p = -(p 2 /h 2 y -g' + 3h 2 (p'/h 2 y + Dh$. (5) 

Figure |l| shows the collapse of a liquid thread of 
propane 6nm in diameter as given by (^J^) with peri- 
odic boundary conditions. The results agree well with 
the computation shown in the supplementary Fig.l of 
JlO|] , and in particular the profile remains close to being 
symmetric. Remembering that P is the velocity of the 
fluid times a typical volume, the multiplicative factor in 
front of the random force in (^|) corresponds to a rela- 
tive increase in noise strength as pinch-off is approached. 
This provides us with a simple physical picture for an 
effective force generated by fluctuations. Namely, a ran- 
dom fluctuation which increases the thread radius also in- 
creases its effective mass, slowing done the motion. Any 
fluctuation towards a smaller neck radius, on the other 
hand, accelerates the motion. On average, the fluctua- 
tions thus drive the thread towards breakup, in fact more 
effectively than surface tension as we will see below. One 
should also note the similarity of (|J) with the equation 
describing directed percolation fll4| , 

FIG. 1. The computed interface profile of a collapsing 
bridge of liquid propane at 185 K. All lengths are nondimen- 
sionalized by the initial bridge radius of 3 nm. The time 
interval between the profiles is 150 ps. 




It is clear that conventional perturbative expansions 
around the deterministic solution are doomed to fail in 
describing this mechanism, since the noise makes a dom- 
inant contribution. In fact, any average taken at a fixed 
time no longer makes sense for this problem, because 
there is a finite probability for pinch-off to have occurred, 
so the original formulation ceases to be valid. Thus a 
valid description has to be conditioned on the event of 
breakup to take place at a fixed time to. It is then 
natural to ask for the most probable sequence of profiles 
that brings one from the initial condition to some final 
profile at time to, corresponding to a "typical" breakup 



event. From a standard path integral description, us- 
ing the Gaussian properties of £, one finds || that the 
probability for such a path is 



W{h(z,t),P(z,t)} ~ exp{- / dtL}, (6) 

Jo 



where the "Lagrangian" is 



L 



2D 2 



J {P + (p 2 /h 2 y + a -3h 2 (p'/h 2 y) 2 _ 

dz - 2 . (7) 



The first equation (^) has no noise in it and has to 
be treated as a constraint. To find the most probable 
path with fixed end-points one derives the Euler-Lagrange 
equation fill for the variation in a = h 2 and P, with the 
constraint treated by adding a Lagrange multiplier d(a + 
P") to L. It is somewhat more convenient to pass to a 
"Hamiltonian" description, introducing P = dL/P as the 
conjugate field, in the literature on critical phenomena 
also known as the "response" field Fr om this one 

directly finds the Hamiltonian equations 



d t a = 
d t P = 

d t a = 

d t P = 



-P" 

-P' 2 /a + 3a(P'/a)' + D 2 aP 
-^-P 2 - PP' 2 /a 2 - 3(PP')'/ 
-2(PP'/a)' - 3((Pa)'/a)' + a 



(8) 
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The contribution G from surface tension was dropped, 
since it will turn out below that it is sub-dominant near 
the pinch point. Introducing the transformation P = 
D 2 P and a = D 2 d the amplitude D can be scaled out of 
the problem. 

"Optimal" paths such as those described by (|J) have 
recently enjoyed some popularity in the statistical me- 
chanics literature Jll|[nJ. However, there are two con- 
ceptual difficulties associated with (||). The first is that 
the equation for P contains a term that corresponds to 
negative diffusion, so it cannot be integrated forward in 
time. In ]lH an ingeneous yet extremely elaborate com- 
putational scheme was developed to deal with this prob- 
lem, based on an initial guess of the complete evolution of 
the profiles. In subsequent iterations, the physical fields 
were always integrated forward in time, the conjugate 
fields backward in time. A second, perhaps more serious 
problem is that one does not know a priori what the final 
profile is supposed to be, so it has to be guessed. Once a 
solution is found, the probability of finding a given final 
profile can be estimated by evaluating the probability of 
the total path. This evidently amounts to a daunting 
mathematical problem for a complicated system like (|^) . 

Both problems can simultaneously be dealt with by 
assuming that the solution is self-similar, as found for 
the deterministic solution [jll 20 1 . This means we assume 
solutions for small \t'\ = \to — t\ to be of the form 
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a(M) = i*'r<m p{z,t) 



l T(£), P(z,t) = \t'\- 1 X (0, 



(9) 



a(z, t) = \t 
i = z/\t'\ 1, \ 

where the exponent a remains to be determined. Plug- 
ging (||) into (||), one finds the similarity equations 

i 
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-2a<t> 2 



A2V 



-2atp 

■4 



r' 



= -r 

= -k-#/^ 4 -3WT/^ 2 



<^ 2 x 



(10) 



x + \x! = -Wir)' - 3((x0 2 )70 2 )' + r", 

where now the prime refers to differentiation with respect 
to £. One notices immediately that (^]) is invariant un- 
der the transformation 6 — > and "0 — > -A 2 ^! so m the 
following it is enough to look for solutions with 6(0) = 1 . 

Physically meaningful solutions of the set (|10|) must 
match onto an outer solution that is slowly varying on the 
local timescale set by \t'\ [0, which means that \t'\ must 
drop out of the similarity description for large arguments 
£ — > ±oo. For the profile <p this means that it must grow 
like 6 ~ £ 2t * at infinity. We are looking for symmetric 
solutions of ([l0|), which were found not to exist for its 
deterministic counterpart |ll|| . Amazingly, we found that 
the subset of solutions of (|10|) with the correct asymptotic 
behavior obey the second order linear equation 



2a i t // 



(11) 



3 T 6 

and the other functions can be written in terms of ip as 
follows: 

</> 2 =0/3, r = -27(V//V) 2 , X = -18(VW, (12) 

which is verified by substitution. Symmetric solutions of 
( |i"T| ) with -0(0) = 3 are also related to the confluent hyper- 
geometric function plj| by ip = 3F(-2a, 1/2, -£ 2 /12). 

FIG. 2. The symmetric similarity profile cf> as given by a 
solution of ©,@ with a = 0.418. 
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This gives a complete description of the physical solu- 
tions of (|l0|), but one still has to find the value of the 
exponent a. To that end we compare the weight (j^) of 
solutions with different a, to find the one for which L is 
minimum. From (0) we find that 



L 



1 



2D 2 



dzaP + Lq, 



where the integration is over the small spatial region 
where the similarity description applies and Lq is essen- 
tially a constant as \t'\ — > 0. 

Keeping the initial height ho of the liquid bridge con- 
stant, we find that the integral J ° L is dominated by 
contributions from \t'\ — > if a < 1/2. Thus in similar- 
ity variables we have up to constants 



hi 



1 



„ dtL ~ 2DH 1 J 2 2a-l/2 



(13) 



The remaining task is to minimize ( |l3| ) as a function of a. 
The decay of the argument of the integral is like £ 4q-4 ; 
so the integral converges for a < 3/4. We conclude that 
( |l3| ) must have a minimum somewhere between a = 1/4 
and 3/4. 

To find it, one has to do the integral numerically for 
general a and finds a m %n ~ 0.418 < 1/2, consistent with 
the assumptions made above. In particular, a is smaller 
than 1, so surface tension becomes asymptotically sub- 
dominant, and hence the singularity is driven by noise 
alone. A more quantitative comparison with numerical 
simulations requires considerably better statistics than 
we are presently able to accumulate. 

We have thus shown that the optimal path method can 
be used to reveal qualitatively new features of nanoscale 
flows, namely a speed-up in pinching and symmetric pro- 
files. It seems that this method is particularly well suited 
to treat problems arising in the new science of nanoscale 
devices. Typically one is interested in the most probable 
or 'typical' behavior that brings one to a certain end, say 
to have a device complete a given motion. I claim that 
the present method is tailored to this situation. 
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